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Abstract. The optical Einstein ring 0047 — 2808 was im- 
aged by the VLT UT1 during the science verification pro- 
gramme. The ring is the image of a high-redshift z = 3.595 
star-forming galaxy, with strong Lya emission at 5589A, 
gravitationally lensed by a massive early-type galaxy at 
z = 0.485. Relative to earlier NTT data the high signal- 
to-noise ratio of the VLT Lya narrow-band image allows 
much improved constraints to be placed on the surface- 
brightness profile of the source and on the mass, leading 
to a measured mass-to-light ratio of M/Ls ~ 13 h for the 
deflector galaxy. We have combined the VLT B-band and 
Lya narrow-band images with a -ftT-band image obtained 
at UKIRT to produce a deep colour image of the system. 
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1. Introduction 

Einstein-ring gravitational lens images should be much 
more common at optical wavelength s than at radio wave- 
lengths (Miralda-Escude and Lehar 1992 ), but so far all 
but one of the known Einstein rings have been discovered 
by radio t echniq ues. The exception is 0047-2808 (War- 
ren et al. 1996a ) where a high-redshift z = 3.595 star- 
forming galaxy, with strong Lya emission at 5589A, is 
lensed by a massive early-type galaxy at z = 0.485. We 
are engaged in a survey to detect similar systems (Warren 
et al. 1996b). The search strategy is to identify anoma- 
lous emission lines (Lya from star-forming galaxies at 
2 < z < 4) in the spectra of a large sample of distant 
early-type galaxies at z ~ 0.4. This has the advantage 



that by the very nature of the identification procedure 
the redshifts of both the source and deflector are obtained 
and so the full lensing geometry is known. In addition, 
because the sources are extended the resulting images, 
rings or arcs, offer the prospect of providing powerful con- 
straints on the m ass distribution in the deflecting galaxies 
(Kochanek 1995| ). Finally, because of the magnification, 
it is possible to study these very fa int sources both spec- 
troscopically (Warren et al. 1998 ) and morphologically, 
resolving angular scales much smaller than is possible for 
unlensed objects. The latter prospects are of particular in- 
terest because the sources are similar to but fainter than 
the population of high-redshift star-forming o bjects iden- 
tified by Steidel and coworkers (Steidel et al. 1996 ), and 



cannot presently be studied in any other way. In this paper 
we present VLT UT1 broad- and narrow-band imaging, 
together with UKIRT if-band imaging of 0047-2808. 

2. Observations and data reduction 

2.1. VLT imaging 

Broad-band B and Lya narrow-band images of 0047-2808 
were obtained on the night of 1998 August 30, with the 
VLT UT1 test camera as part of the Science Verification 
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* Based on observations obtained at the European Southern 
Observatory, Paranal, Chile, and the United Kingdom Infrared 
Telescope, Hawaii 



programme (Leibundgut et al. |1998|). The CCD pixel size 
is 0'.'045. However the CCD was binned 2 x 2, so the pixel 
size in all the frames was 0'.'09. The narrow-band filter has 
a central wavelength 5589A and width 20A FWHM. Inte- 
gration times were 3 x 300 sec (B) and 3 x 1200 sec (Lya). 
The seeing was 0.6—0.7". Procedures followed for bias sub- 
traction and flatfielding were mostly standard. However 
the flatfielded narrow-band images required a correction 
for large-scale gradients. This was achieved by firstly com- 
bining the deregistered frames, clipping out objects. The 
resulting frame was smoothed, and normalised, and the 
flatfielded frames were divided by this correction frame. 
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Fig. 1. True-colour image of the field of the Einstein-ring gravitational lens 0047-2808, showing a region 68" x 81". 
North is up and East to the left. The image is a rgb combination of the VLT B and narrow-band images and the 
UKIRT K image. Because of the strong Lya emission line in the narrow-band the ring (z = 3.595) appears as green. 
The lensing galaxy (z = 0.485) is red because the galaxy is bright in the K band but faint in the B and narrow bands. 
There is a second massive elliptical galaxy in the field, visible as the brightest red image in the SW quadrant, which 
has a redshift z = 0.488. 
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2.2. UKIRT imaging 

Broad-band K images of 0047-2808 were obtained with 
the UKIRT IRCAM3 instrument on the nights of 1997 
September 12 and 13. The pixel size was 0'.'286. The fi- 
nal image is a mosaic from two positions, one centred on 
0047-2808, total integration time 15 120 sec, and another 
at a position 58" to the SSW centred on a second distant 
elliptical, total integration time 4725 sec. At each position 
several sequences of 9-point dithers were summed. The 
seeing averaged 0'.'7. The data were flat -fielded using a se- 
quence of twilight sky exposures, and then an appropriate 
sky frame, formed from a running median filter through 
the stack of images, was subtracted from each data frame, 
and the resulting frames registered and summed. 

2.3. Results 

Fig. 1 shows the rgb colour image resulting from combining 
the B (=&), Lya (—g), and K (—r) images. The ring stands 
out strongly in green because of the strong Lya line in 
the narrow-band filter, while the lensing galaxy is very 
red and is visible inside the ring. A minimum % 2 fit of a 
de Vaucouleurs model for the light profile of the lensing 
galaxy in the fsT-band image was computed by convolving 
two-dimensional r 1 / 4 profiles with the psf, measured from 
a star in the frame. (The if -band image is best for fitting 
the galaxy profile because the ring is not detected at this 
wavelength, and the contrast between the galaxy and the 
ring is maximised.) The model was then convolved with 
the Lya-band psf, scaled to the central counts in the Lya 
image, and subtracted. The resulting image of the ring, 
rcbinncd to a pixel size of 0'.'18, is shown in the top left- 
hand panel of Fig. 2. 

3. Gravitational lens model 

Comp ared with the original NTT image (Warren et al. 
1996a ) the new VLT image has much higher signal-to- 
noise ratio. The counter image predicted by our original 
model, but not convincingly detected in the NTT image, 
is now clearly seen. The same modelling procedure as de- 
scribed in Warren et al. (1996a), where the projected sur- 
face mass density was assumed to follow the (intrinsic) 
de Vaucouleurs profile, now measured from the if -band 
image (as described above), has been applied. The single 
free parameter is the global mass-to-light ratio (M/L). 
Utilising the computational technique for arbitrary 



where V = 2d/dz. The deflection potential, is related 
to the projected surface mass density in the lens, E, by 



lenses with elliptical symmetry (Schramm 1994) the M/L 
ratio in the model was adjusted to produce the most com- 
pact configuration for the unlensed image in the source 
plane. Here we briefly review the key steps in the proce- 
dure; a position in the image plane, represented by the 
complex coordinate z, is mapped onto the source plane 
position to, according to 



V 2 $(z,z) = 2E(z,z)/E crit 



(2) 



where E cr j t is the critical surface mass density to gravita- 
tional lensing and is given by 
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w(z,z) — z - V$(z,z) 



(1) 



Here, Dy are the angular diameter distances betwe en th e 
source (s), lens (1) and observer (o) (Schneider et al. 1992 ). 

A 61x61 pixel (5.5 x 5.5") region of the VLT narrow- 
band image, centred on the galaxy, was used for the com- 
putation. Assuming a fiducial value of M/L, the mapping 
given in Equation [j] gives the coordinates in the source 
plane of any image-plane coordinate. The M/L was ad- 
justed, focusing the emission over the source plane into a 
small region. This determines the centroid of the source. 
The source was then modelled as a Gaussian profile. The 
structure in the ring (i.e. the angular extent of the gaps, 
size of the counterimage) is dictated by the angular extent 
of the source. A source of FWHM of 0"2 when reimaged 
by the lensing potential was found to reproduce the struc- 
ture in the ring well. To reimage the source each pixel was 
sub-pixelated into a 10x10 grid; these grid points were 
mapped to the source plane to measure the surface bright- 
ness at each grid point. Mapping of the surface brightness 
in this way is accurate provided the grid spacing mapped 
to the source plane is substantially smaller than the scale 
over which the surface brightness of the source varies. 

Having fixed the source position and profile the lens 
M/L was then finely readjusted to provide the best fit, in 
terms of \ 2 , of the model of the ring to the data. The re- 
sults of this procedure are presented in Fig. 2. The upper 
left-hand panel presents the VLT image of the ring after 
subtraction of the model for the surface brightness distri- 
bution of the foreground galaxy. Below, the model source 
is shown on the same scale together with the caustic lines 
defined by the gravitational lens model; the caustics de- 
lineate the regions of multiple imaging over the source 
plane. The resultant image for this source configuration 
is presented in the lower right-hand panel. In the upper 
right hand panel this image has been convolved with a 
Gaussian seeing profile to approximate the observing con- 
ditions. There is good correspondence between the struc- 
ture in the model and the observed structure of the ring. 

The measured angular radius of the ring in the VLT 
image is 1'.'08 ± 0.03. This more accurate value is smaller 
than the value measured from the NTT image (1" 35 ±0.1) 
by 20% and this significantly lowers the mass estimate. 
Part of the discrepancy between the two measurements 
is due to the fact that the ring is elliptical in shape and 
that the counterimage (invisible in the old data) lies on 
the minor axis i.e. the old value for the radius was mea- 
sured along the major axis of the ellipse. The two measure- 
ments are consistent therefore. The computed mass within 
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Fig. 2. False-colour images displaying the results of modelling the mass distribution in the lens. The size of each square 
image is 5.5" on a side. North is up and East to the left. The top LH image shows the VLT Lya image after subtraction 
of the image of the lensing galaxy, and rebinned to a pixel size of 0'.'18. The lower LH image shows the model of the 
source, and the 3-image and 5-image caustics. The lower RH image shows the image that would be observed at infinite 
spatial resolution, and the upper RH image is the lower RH image convolved with the psf to emulate the observing 
conditions. The upper RH image therefore models the upper LH image with zero noise. The good correspondence of 
the two upper images gives a measure of the accuracy of the mass model. 
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the Einstein radius is 1.73(1.95) x lO n /j _1 M for q Q = 
0.5(0.1) The uncertainty in the mass estimate within 
the Einstein radius is dominated by the uncertainty in the 
radius rather than the form of the mass profile. Changing 
the radius by C'06 (i.e. 2cr) changes the computed mass 
by 10%. The M/L ratio for the model, corrected for lumi- 
nosity evolution (Paper I), is M/Lb(o) = 14.2 ft,(12.1 h). 
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1 h is the Hubble constant in units of 100 km s *Mpc 
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